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Abstract 

I give a pedagogical account of Shor's nine-bit code for correcting arbitrary errors on single 
| qubits, and I review work that determines when it is possible to maintain quantum coherence 

by reversing the deleterious effects of open-system quantum dynamics. The review provides an 
opportunity to introduce an efficient formalism for handling superoperators. I present and prove 
some bounds on entanglement fidelity, which might prove useful in analyses of approximate error 
correction. 



1 Introduction 

A bit is the fundamental unit of information, represented by a choice between two alternatives, 
conventionally labeled and 1. In the real world the abstract notion of a bit must be realized as a 
physical system. A classical bit, the unit of classical information processing, can be thought of as 
a two-state classical system. The two states, perhaps or 1 printed on a page or the two positions 
of a particle in a double-well system, are distinguishable, and because they are distinguishable, 
^ ■ they can be copied. A quantum bit or qubit, the unit of quantum information processing, is a two- 

state quantum system. The two basis states, labeled |0) and |l), are distinguishable and can be 
copied, just like the states of a classical bit. The difference arises from the superposition principle: 
the general state of a qubit is an arbitrary linear combination of |0) and |l). The many possible 
superposition states available to a qubit give it more information-processing power than a classical 
bit, even though the general superposition states cannot be distinguished reliably and cannot be 
copied — the no-cloning theorem forbids |l[E|. The enhanced information-processing power of qubits 
can be harnessed to a variety of information-processing tasks including computation (for a review 
and extensive list of references, see ||). 

The price for the power of quantum information is eternal vigilance in maintaining quantum 
coherence, for the enhanced information-processing power comes from the ability to manipulate 
superposition states. Coupling to the environment, with the accompanying noise and decoherence, 
tends to destroy superpositions. Indeed, what we call a classical bit is just a qubit whose coupling 
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to the environment keeps it from occupying superpositions of two orthogonal states. Decoherence 
strips the kets off |0) and |l), leaving the and 1 of a classical bit, which cannot be superposed. 

Quantum information processing got a tremendous boost with Shor's surprising discovery ||| 
that quantum information stored in superpositions can be protected against decoherence. Shor's 
announcement of a nine-bit quantum code, followed shortly by the discovery of five-bit || and seven- 
bit codes ||, ignited an explosion of activity on quantum error correction ||, || [H], |ll|, 12] (see 
[||] for further references), which led to the demonstration that a quantum computer can perform 
arbitrarily complicated computations provided that the error per operation can be reduced below 
an error threshold gj| 0, [D| . 

In this article I give in Sec. ^ a pedagogical account of the simplest error-correction scheme, 
Shor's nine-bit code ||, with the aim of illustrating the essential ideas of quantum error correc- 
tion. In Sec. |3| I review work that determines when it is possible to maintain quantum coherence 
by reversing the deleterious effects of open-system quantum dynamics, using the review as an op- 
portunity to introduce an efficient formalism for manipulating superoperators. In Sec. || I first 
review the information-theoretic description of error correction and then present and prove some 
bounds on entanglement fidelity, which might turn out to be useful in considering approximate 
error correction. 



2 Quantum Error Correction: Shor's Nine-Bit Code 

Quantum error correction is closely related to classical error correction. Much of the theory of 
quantum error correction [7L |8|, ^] comes from the theory of classical linear codes, in which a classical 
code word of length k, i.e., a string of k O's and l's, becomes a vector in a /c-dimensional vector 
space over the field consisting of and 1. I do not discuss this theory here, but rather consider 
a particular quantum code, Shor's nine-bit code, to illustrate the essential ideas of quantum error 
correction. 

To get started, though, let's first consider correcting errors on a classical bit. The only type of 
error is a bit flip, in which and 1 are interchanged. If the errors are rare, they can be corrected 
using redundancy, because one only needs to correct single-bit errors. Suppose, for example, that 
a is encoded as the three-bit string 000, and a 1 as the three-bit string 111. Then an error on a 
single bit can always be detected and corrected. The situation is summarized here: 

no error 000 111 , 

flip 1st bit 100 011 , , . 

flip 2nd bit 010 101 , ^ ' 

flip 3rd bit 001 110. 

An error can be detected by polling the three bits. The minority bit suffered the flip, which can 
be corrected by flipping the minority bit again. If the probability for a bit flip on a single qubit 
is p <C 1, then this scheme reduces the probability of error from p to 3p 2 , a winning proposition if 

The key elements in this scheme, redundancy and polling, cannot be translated to quantum 
error correction: making redundant copies of a qubit is ruled out by the no-cloning theorem, which 
forbids making copies of arbitrary superposition states; polling requires ascertaining the state of 
each qubit, which destroys the quantum coherence one is aiming to protect. On closer examination, 
though, the situation is more promising. The eight strings that result from no error and from the 
three bit flips are all distinct. This distinguishability opens up the possibility of a quantum code 
in which |0) is encoded as a state of three qubits, the "logical zero" state \0)l = |000), and |l) 
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is encoded as the "logical one" state = |m)- The logical zero and logical one span a two- 
dimensional subspace, the "code subspace" to be used for quantum information processing. Notice 
now what happens to an arbitrary superposition of the logical one and logical zero states under 
single-qubit bit flips: 



no error 1(8)1(8)1 a|000) + j3\ 111) 

flip 1st qubit o\ (8 1 <8 1 ajlOO) + /5|011> 

flip 2nd qubit 1 <8 ai (g> 1 a|010) + (3\ 101) 

flip 3rd qubit 1 <8 1 <8 <Ji ajooi) + /3| 110) 



(2) 



Here the middle column writes the error in terms of the unit operator 1 and the bit-flip Pauli 
operator o\ = |0)(l| + |l)(0| for the appropriate qubit. 

The no-error operator and the three bit-flip errors map the code subspace unitarily to orthog- 
onal two-dimensional subspaces within the eight-dimensional Hilbert space of a qubit triplet. The 
single-qubit errors can be detected and distinguished by determining in which two-dimensional 
subspace the system lies, without in any way disturbing the superposition, and the error can be 
corrected by mapping the error subspaces unitarily back to the code subspace. This is the essence 
of quantum error correction: find a code subspace such that the high probability errors map uni- 
tarily to orthogonal subspaces; then the errors can be detected and corrected without destroying 
quantum coherence. The superposition states in the code subspace are entangled states of the three 
qubits. The redundancy used by a classical code, which cannot be translated to quantum coding, 
is replaced by entanglement in a quantum code. 

Having gotten this far, however, we now realize that the task is tougher than classical error 
correction, because there are quantum errors that have no classical counterpart. Specifically, there 
are "phase flips," described by the Pauli operator 0-3 = |0)(0| — |l)(l|, and errors described by the 
Pauli operator — 102 = C103 = — |0)(l| + |l)(0|, which is a phase flip followed by a bit flip. If we 
can correct these errors, we can correct all single-qubit errors, because all error operators can be 
written as a linear combination of the unit operator (no error) and the three Pauli operators. 

Let's concentrate first on the phase-flip errors, hoping that the combined phase-bit flips take 
care of themselves. The first thing to notice is that o\ and 0-3 switch roles in the transformed basis 
defined by 

|±) = -^(|0)±|1» : (3) 
<t 3 becomes a bit flip, and o\ becomes a phase flip, 

<7 3 |±) = |=F) , <7!|±) = ±|±). (4) 

Thus we can correct single-qubit phase flips by using a quantum code whose logical basis states are 

|+++) and I ), but this comes at the expense of being unable to correct the original bit-flip 

errors. 

Here entanglement comes to the rescue again. The entire code subspace spanned by |000) and 
I 111) is protected against single-qubit bit flips, so we can use any orthogonal basis in the code 
subspace as the logical zero and one. In particular, we can use "up" and "down" states 



|T>^(|ooo> + |iii», 

|i> = ^(|000>-|iii», 



(5) 

which are defined in analogy to the way the states |±) are related to |0) and |l) for a single qubit. 



3 



The effect of single-qubit bit flips on these states is summarized here: 
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The eight types of up and down states make up an orthogonal basis for the three-qubit Hilbert 
space. 

Notice now that the bare up and down states are flipped by all three single-qubit phase flips: 
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This suggests correcting both bit- and phase-flip errors by again tripling the number of qubits and 
using the following nine-qubit states as the logical zero and one states: 







ITTT) 



HID 



(8) 



This logical zero and one constitute Shor's nine-bit code. 

To complete the discussion of the nine-bit code, I display in Table p] how an arbitrary state 
in the code subspace is affected by all 27 single-qubit errors. Notice first that in accordance with 
the discussion above, the three phase-flip errors on each qubit triplet have exactly the same effect. 
Thus the nine-bit code is a degenerate code ||, one in which errors that are independent on the 
entire Hilbert space become the same on the code subspace. The degeneracy leaves a situation 
where we must consider the no-error case plus 21 independent errors. Each error maps the code 
subspace unitarily to a two-dimensional subspace. Moreover, by examining the table, one sees that 
the code subspace and the 21 error subspaces are mutually orthogonal. The 22 states that arise 
from |0) l are orthogonal to the 22 states that arise from |1)l : the former all have up states for two 
of the triplets, whereas the latter all have down states for two of triplets, so they disagree on up 
versus down in at least one position. The states that arise from \0)l (|1)l) are mutually orthogonal 
because they are the 22 states that come from putting two of the triplets in the up (down) state, 
with the third triplet cycling through the eight up- and down-type states that span the triplet 
Hilbert space. Single-qubit errors can be corrected by determining in which of the orthogonal 
subspaces the nine qubits lie and then mapping that subspace unitarily back to the code subspace. 

As noted above, error correction works for single-qubit errors that are not described by Pauli 
operators, because any error can be written as a linear combination of the unit operator and 
the Pauli operators. Consider, for example, decay from |0) to |l), with essentially instantaneous 
phase decoherence between |0) and |l), a situation described by three error operators: A% = 
VT=t|0)(0| = y/T=i(l + o 3 )/2, A 2 = V7|l)(0| = ^(01 - ia 2 )/2, and A 3 = |1)<1| = (1 - <r 3 )A 
7 being the probability of decay. If one of the nine qubits suffers a decay, a measurement of the 
error subspace reveals either no error or one of the Pauli errors on that qubit. All of these being 
correctable, it doesn't matter which is the result of the measurement. 
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The nine-bit code wastes Hilbert space, for it uses only 44 of the 2 9 = 512 dimensions in the 
nine-qubit Hilbert space. If one wants to correct r errors per qubit, using a code with iV qubits, 
then to accommodate the code subspace and the rN errors, one needs 2(1 + rN) dimensions. This 
leads to the quantum Hamming bound M: 

2 N-i _ l 

2(l + rN)<2 N r< — . (9) 

Three qubits permit correction of one error, the situation we started with in this section. Five 
qubits have the potential for correcting the three Pauli errors, a potential realized in an astonishing 
five-qubit quantum code [||] . 



3 Reversal of Open- System Dynamics 

In this section I review the description of open-system dynamics in terms of quantum operations and 
the question of when quantum coherence can be maintained by reversing open-system dynamics. 
The review serves to introduce a formalism for handling quantum operations, which provides insight 
into their structure. 

Throughout this section and the next I use a set of conventions introduced by Schumacher 
[16]. The primary quantum system is denoted by Q; it is assumed to be finite-dimensional, with a 



Hilbert space TLq of dimension D. The primary system interacts with an environment E, and to 
deal with purifications of Q states, there can be a additional, passive reference system R. Where 
it is necessary to avoid confusion, superscripts R, Q, and E are used to distinguish states and 
operators of these systems. Initial states are unprimed, and states after the dynamics are denoted 
by a prime. 



3.1 Open-system dynamics and quantum operations 

Consider a primary system Q, initially in state p, which is brought into contact with an environment 
E, initially in state p E = Aj|0j)(0j|, where the states \4>i) are the eigenstates of p E . The 
two systems interact for a time, the interaction described by a unitary operator U, and then the 
environment is observed to be in a subspace spanned by orthogonal states \gj), corresponding to a 
projection operator P E = J2k \9k}(dk\- The unnormalized state of the system after the observation 
is given by a partial trace over the environment, 

tv E (P E U{p ® p E )tf) = A{p) , (10) 

where A is a linear map on system density operators. Inserting the forms of the projector P E and 
the environment state p E leads to 



A(p) = VXii3k\U\4>i)p{<t>i\U^\9k)VXi = £ A*p4 • (11) 

k,l a 

The system operators 

A a = A kl = ^\ l {g k \U\<Pi), (12) 

where the Greek index a is as an abbreviation for k and I, provide an operator decomposition of the 
map A and thus are called decomposition operators. The normalized post-dynamics system state 
is 

P = W At \\ ' 13 ) 
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where one should notice that 



tr(P E U(p ® p E )tf) = tx(A(p)) = tr( P J2aU«) (14) 

is the probability for the environment to be found in the specified subspace. 

A primary quantum system that is exposed to an initially uncorrelated environment always 
has dynamics described by a map like A. This includes both nonselective dynamics, where no 
observation is made on the environment (P E = 1 E ), and selective dynamics, where the system 
state is conditioned on the result of a measurement on the environment (P E ^ 1 E ). It is important 
to identify the mathematical conditions that characterize a suitable map A. We can immediately 
identify three conditions that A must satisfy. 

Condition 1. A is a linear map on operators; i.e., it is a superoperator. 

Condition 2'. A maps positive operators to positive operators. (A positive operator G is one 
such that (tp\G\tp) > for all vectors \tp); density operators are positive operators.) 

Condition 3. A is trace- decreasing, i.e., tr(A(p)) < 1 for all density operators p. This condition, 
which follows immediately from Eq. (|l~4l), can be expressed as an operator inequality, 

J2 A Ua<L (15) 

a 

The trace is preserved for nonselective dynamics, but generally decreases for selective dynamics. 

It might be thought that the above three conditions are sufficient to characterize A, but it 
turns out that Condition 2' must be strengthened. The more restrictive condition can be motivated 
physically. Suppose that R is a reference system that, though it does not take part in the dynamics, 
cannot be neglected because the initial state p of Q is the partial trace over R of a joint state p R ® . 
We certainly want the map I R <S> A, where I R is the unit superoperator on R, to take p R ® to a 
positive operator, which can be normalized to be a density operator. This condition holds trivially 
for a map of the form (|Tl| ) , 



{1 R ® A)(p RQ ) = J2(1 R ® A a )p RQ (l R <g> At) > , (16) 

a 

so we replace Condition 2' with a stronger condition. 

Condition 2. A is completely positive; i.e., (T R (g) A)(p R ®) > for all joint density operators 
pRQ Q £ q anc j arD itrary reference systems R. 

A map on operators that satisfies Conditions 1-3 is called a quantum operation. The description 
of open-system dynamics in terms of quantum operations was pioneered by Hellwig and Kraus 



[17, 18, 19 1 . Thus far our discussion has established that any open-system dynamics of the sort 



introduced above is described by a quantum operation. We need two further properties: first, 



that any quantum operation has an operator decomposition, as in Eq. (11), and second, that any 
quantum operation can be realized by a unitary coupling to an environment. Once the former 
result is established, the latter is easy. Any linear, trace-decreasing map that has an operator 
decomposition can be realized by a unitary coupling to an initially pure-state environment. One 
partially defines a joint operator U on RQ as in Eq. ( |i"2| ) and uses the trace-decreasing condition ( |l5"| ) 
to show that the definition can be extended so that U is unitary. We turn now to the connection 
between complete positivity and the existence of an operator decomposition, which requires us to 
step back and consider the properties of superoperators. 
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3.2 Superoperators and complete positivity 

The space of linear operators acting on Hq is a Z) 2 -dimensional complex vector space C(Hq). Let 
us introduce operator "kets" \A) = A and "bras" (A\ = A< , distinguished from vector kets and bras 
by the use of smooth brackets. Then the natural inner product on C(TCq), the trace-norm inner 
product, can be written as (-A|2?) = tr(^4'i?). An orthonormal basis \ej) induces an orthonormal 
operator basis | e^) (e^ | = Tj k = r a , where the Greek index is again an abbreviation for two Roman 
indices. Not all orthonormal operator bases are of this outer-product form. 

The space of superoperators on Q, i.e., linear maps on operators, is a D 4 -dimensional complex 
vector space C(C(TLq)). Any superoperator S is specified by its "matrix elements" 

<%, m jfc = (ei\S{\ej)(e k \)\e m ) , (17) 

for the superoperator can be written in terms of its matrix elements as 

S= ^ Si j>mk \ei)(e j \(g)\e k }(e m \=^S a pT a (g)Tj j = ^2S a p\T a )(Tf } \ . (18) 

lj,mk a,/3 a,f3 

The ordinary action of S on an operator A, used above to generate the matrix elements, is obtained 
by dropping an operator A into the center of the representation of S, in place of the tensor-product 
sign, 

S{A) = Y,S aP T a ATl . (19) 

There is clearly another way that S can act on A, the left-right action, 

S\A) = J2^\ra)(rp\A), (20) 

in terms of which the matrix elements are 

S af 3 = (T a \S\Tp) = (Tij\S\T rnk ) = (ei\S(\ej)(e k \)\e m ) . (21) 

This expression provides the fundamental connection between the two actions of a superoperator. 

With respect to the left-right action, a superoperator works just like an operator. Multiplication 
of superoperators T and S is given by 

TS = ]T T ai S lP \T a ){T P \ , (22) 

and the adjoint is defined by 

(A\tf\B) = (B\S\AT ^ S^ = J2^ a \r a )(rp\. (23) 

With respect to the ordinary action, superoperator multiplication, denoted as a composition To S, 
is given by 

T o S = ]T T l5 S aP r 7 r a t\t\ . (24) 

The adjoint with respect to the ordinary action, denoted by S x , is defined by 

tr([5 x (B)]U)=tr(Bt 5 (A)) ^ S» = ® r . (25) 

a,/3 
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To deal with complete positivity, we need to introduce a reference system R, which we choose 
now and henceforth to have the same dimension as Q, and we need to have a way of turning oper- 
ators (superoperators) on Q into vectors (operators) on RQ. To do so, introduce the unnormalized 
maximally entangled state 

l*) = El/i>®l e i> = El/i> e ;>> ( 26 ) 

j i 

where the vectors |/,-) comprise an orthonormal basis for R. The VEC of an operator A on Q is 
the vector 

\$ A ) = 1 R ® A\V) = ]T \fj) ® A\ ej ) . (27) 

3 

The VEC map is a one-to-one, linear map from C(TCq) to TCrq; we recover A from \$>a) y i a 

{fj,e k \$ A ) = (e k \A\ ej ) . (28) 
It is easy to see that VEC preserves inner products, 

($a\$b) = tr(AtB) = (A\B) . (29) 
One further aspect of VEC deserves mention. Given a density operator p for Q, applying VEC to 

|<^> = i R ^VpI*) = E l/i) ® V^l e i) , (3°) 

3 

generates a purification of p, i.e., 

tr H (|^)(^|)=p. (31) 

The analogous OP map is a one-to-one, linear map from C{C{TLq)) to C(TLrq). It takes a 
superoperator 5 on Q to the operator 

(Z*®S)(|*)(v|/|) = ^2 \fj)(fk\ ®S(\ ej )(e k \) , (32) 
3,k 

and we recover 5 via 

(33) 

Looking at OP in a slightly different way, 

(i R ®s)mm=j2saM*T a )(* Tf3 \ , (34) 

we find that 

(<S> A \(l R ®S)(\*)m\<S>B) = (A\S\B) . (35) 

Thus the OP of a superoperator S operates in the same way as the left-right action of S. 

We are ready now to return to complete positivity. Recall that we are trying to show that 
any completely positive superoperator A has an operator decomposition. The key point is that 
complete positivity requires that the OP of ^4, i.e., (I (g) ,A)(|\E'}(\E r |), be a positive operator, but 
Eq. (|35| ) now shows this to be equivalent to the requirement that A be positive relative to its 
left-right action, which I write as A > 0. Any such positive superoperator has many operator 
decompositions, including its orthogonal decomposition. Moreover, we get for free, by using the 
decomposition theorem for positive operators pOfl, the following result, originally due to Choi plf: 
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two decompositions A a and B a give rise to the same completely positive superoperator if and only 
if they are related by a unitary matrix Vg a , i.e., 

B/3 = J2 V ^A a (36) 

a 

(if one decomposition has a smaller number of operators, it is extended by appending zero opera- 
tors). 

What we have shown in this subsection is that a map is completely positive if and only if it is a 
positive superoperator relative to the left-right action. For a quantum operation we must add the 
trace-decreasing condition (|i~5|), which now can be put in the compact form 

A*(l)=J2AlA a <l, (37) 

a 

with equality if and only if the operation is trace preserving. Thus we can now characterize a 
quantum operation as a superoperator that is positive relative to the left-right action (.4 > 0, 
complete positivity) and that satisfies -A x (l) < 1 (trace-decreasing). 

We should introduce one more ingredient before leaving our discussion of superoperators. Sup- 
pose the initial state p of Q is VEC'ed to an initial joint state |$ ^p-), as in Eq. (|30|), The joint state 
of RQ after the dynamics described by T R (g> A is given by 

9 ~ V Z ® tr(A(p)) ■ (38) 



Referring to Eq. (|35|) , we see that the superoperator 

_ Ao^p®^p 
Ap ~ tr{A{ P )) 1 ' 

is equivalent to the joint density operator p R ®' ; i.e., pPQ' operates the same as the left-right action 
of A p . 

Straightforward consequences of the definition of A p are that 

A ^ = TO = ^' <*» 

P V ^ tr(.4(/>)) " tr(,4(p)) V ; 

If .4 is trace preserving, then the density operator a = p. The physical significance of a can be 
ferreted out with a bit more work. First we note that 

(e k \a\ ej ) = tr(^ p x (l)| ej )(e fc |) = tv(A p (\ ej ) (e k \)) . (42) 

Writing the joint state of RQ after the dynamics as 

ft = (i R ®A p )mm = J2\fj)(fk\®A p (\ ej )(e k \) , (43) 

j,k 

we find that the state of R after the dynamics, 

p Rl = tv Q (p R Q') = £ \fj)(fk\(ek\"\ej) , (44) 

is the "transpose" of a with respect to the bases \ej) and \fj). 
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3.3 Reversal of quantum operations 

We are now ready to formulate the problem of reversing open-system dynamics — i.e., correcting 
errors due to coupling to an environment. In this subsection I follow closely the formulation found 
in p3| . The mathematical statement of the problem is the reversibility of a quantum operation 
A on a "code subspace" C of the system Hilbert space TCq. The reversal must be accomplished 
by a physical process, so it, too, is described by an operation, 1Z. We want the reversal definitely 



to occur, so we require that 7Z be a trace-preserving operation. Thus we say that |1C], 22, 23] a 
quantum operation A is reversible on the code subspace C if there exists a trace-preserving reversal 
operation 1Z, acting on the total state space of Q, such that 

, K o A(p) 

np) = l4A(p)j = P (45) 

for all p whose support is confined to C. An immediate consequence of the linearity of 1Z o A is 
that 

tr(A{p)) = constant = p 2 (46) 

for all p whose support is confined to C p2| ], where p is a real constant satisfying < p < 1.. This 
allows us to rewrite the reversibility condition (|45|). To do so, we introduce the restriction of A to 
C, 

Ac = AoP c ®P c , (47) 
where Pc is the projector onto C. Then the reversibility condition becomes 

K o Ac = p 2 Pc ® Pc ■ (48) 

It is not hard to show (for proofs see fi~l"l , |23|) that a quantum operation A, with decomposition 
operators A a , is reversible on C if and only if there exists a positive matrix m a p, having unit trace, 
such that 

P c A\A a P c = p 2 m a(3 P c . (49) 

It is instructive to rewrite this condition as (ej\A^A a \ek) = p 2 m a p5jk, where the vectors \ej) make 
up an orthonormal basis on C. We can think of each decomposition operator A a as an "error 
operator." Though each error operator acts like a multiple of a unitary operator within C, different 
error operators are not required to map to orthogonal subspaces. How do we square this with the 
discussion in Sec. |2|? 

The puzzle is resolved by realizing that the decomposition of A is not unique. We need to 
choose the decomposition operators to represent independent, indeed orthogonal errors within C. 
To do so, take any density operator p whose support is the entirety of C, and use a unitary matrix 
transformation of the decomposition operators A a to diagonalize the matrix p 2 m a p = t^{A a pJ^ ^ 



In the new decomposition, called a canonical decomposition [|23| , the decomposition operators A a 
satisfy 

' ' 1 ~ - 1 "H> 

where the eigenvalues X a satisfy 1 = J2 a X a . In terms of this canonical decomposition, the reversal 
condition (H9f) becomes 



(Ai3^/p\A a ^/p) = tr(A a pAi) = p 2 m Ql3 = p 2 \ Q 6 a[ 3 , (50) 



P c A\A a P c = p 2 X a 5 a(3 P c ■ (51) 
This condition has a ready interpretation. When a = j3, the operator polar-decomposition 



theorem 24] implies that there exists a unitary operator U a such that 



i a P C = U a \JP C AiA a P C = pVXaU a Pc = P^\ a P a U a . (52) 
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In the last equality we introduce the projector P a = UaPcU^ onto the subspace that is the unitary 
image of C under U a . In terms of the canonical decomposition the restriction of A to C becomes 

Ac = E A « P C ® P C^i = E X * U * P C ® P £4 = ^ E AaPa^a ® ■ (53) 

a a a 

When one or more of the eigenvalues A Q is zero, we are dealing with a degenerate code ||, 23 1: as 
far as the action of A within C is concerned, the X a = operators are irrelevant. We discard those 
operators henceforth, remembering that this is legitimate so long as we restrict attention to the 
code subspace. The irrelevant decomposition operators discarded, the content of the a ^ (3 terms 



in Eq. (49) is that the image subspaces P a are orthogonal, i.e., 

P*Pp = SapPp ■ (54) 

Thus in the canonical decomposition, the error operators A a act like multiples of unitary operators 
on C and map C to orthogonal subspaces. 

It is important for the considerations in Sec. to note that when we go to a canonical decom- 



position, we are diagonalizing the completely positive superoperator (39|), 



A p = E^WP®VMa = ^KU ay /p®y/pUl =Y, X <*VP^ U a®Ul^ , (55) 

a a a 

for the operators A a ^J~p are orthogonal according to Eq. (j5C|). In the last equality of Eq. (p5|) , 
we introduce the orthogonal density operators p a = U a pU^ that p is mapped to by the unitaries 
U a . Notice that the eigenvalues and normalized eigenoperators of A p (relative to the left-right 
action) are X a and U ay /p = y/p^U a . The equivalence between A p and p R ® means that OP'ing the 
operators U ay fp generates the eigenvectors of p R ®' . 

For our purposes the only relevant part of the reversal operation is its restriction to the subspace 
N that is the direct sum of the unitary images of C under the unitaries U a . This restriction, which 



is unique |23], is given by 

n N = Y,UtP a ®P*U a . (56) 

a. 

It is easy to verify that this 1Z n reverses Ac, and it also satisfies TZ^(1) = ^2 a P a = Pn, the 
appropriate trace-preserving condition for the restriction. It is interesting to note that since 

p' = E X aU a p t4 = E X aPa (57) 
a a 

can be written in terms of an ensemble of orthogonal density operators p a , we can give a very 
compact equation for the relevant part of the reversal operation: 

K p . = A* . (58) 

4 Entropy Exchange and Entanglement Fidelity 
4.1 Information-theoretic formulation of reversibility 



Our starting point in this section is Eq. (|57|), which gives p' in terms of an ensemble of orthogonal 
density operators p a . This allows to conclude that 

S(p') = 5(p)-EA a logA Q , (59) 
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where S(p) = —ti(plogp) is the von Neumann entropy of p. Since the quantities X a are the eigen- 
values of A p , the second term on the right is the entropy of the completely positive superoperator 
A p or, equivalently, the entropy of p R ® . Schumacher [16| introduced this entropy and dubbed it 
the entropy exchange 

S e (p, A) = S(p RQ ') = S(A P ) = - £ Aa log A Q . (60) 

a 

It is useful at this point to introduce the superoperator trace 

Tr(S) = 5>a| S\r a ) = tr((l R ® = tr(<S(l)) , (61) 



where the operators r a make up an orthonormal operator basis and where we use Eqs. ( p2| ) and 
to reduce the definition to an operator trace. Not surprisingly, A p has unit trace: 

Tr(A p ) = tv(A p (l)) = tr(p') = 1 . (62) 



The entropy exchange (60) can now be written as 

S(A P ) = -Tr(A p logA p ) . (63) 
Equations (|46] ) and J59| ) are consequences of the reversibility of A on C. It turns out that they 



are also sufficient to ensure reversibility [23|. A quantum operation A is reversible on the code 
subspace C if and only if the following two conditions are satisfied: 
Condition 1. 

HA(p)) = p 2 (64) 

for all p whose support is confined to C, where p is a real constant satisfying < p < 1; 
Condition 2. 

S(p) = S(p') - S(A P ) (65) 

for any one p whose support is the entirety of C ( and then for all p whose support is confined to 
C). A proof of this theorem is given in 23]. Here I give a very simple proof of sufficiency, necessity 
already having been demonstrated. 

We start with the orthogonal decomposition of 

A P = Y, X *T a ® T-t , (66) 
a 

where the eigenoperators are orthonormal, i.e., (r Q |r^) = 5 a p. The polar-decomposition theorem 
[24] guarantees that there exists a unitary operator U a such that 

T a = Ua^/da = yfp^U a , (67) 

where a a = T^T a and p a = T a T^ are normalized density operators that are unitarily equivalent, i.e., 
Pa = UaOaU^. The unitary equivalence implies that S(p a ) = S(a a ). We have from Eqs. ( |40| ) and 
@ that 



p' = A p (l) = ^XaPa, 

a 

a = A x p {l) = Y. X « a <*- ( 69 ) 
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These ensembles for p' and a give several inequalities |25|, p6[ : 

S(A P ) > S(p')-J2^S(p a )>0, (70) 

a 

S(A P ) > S(a)-J2^S(a a ) >0. (71) 

Equality holds on the left if and only if the density operators p a (a a ) are orthogonal, whereas 
equality holds on the right if and only if the density operators p a = p' (a a = a) for all a. 

We have not yet used the two Conditions. If write Condition 1 as p? = tr(„4(p)) = tr(_4 x (l)p) 
for all p whose support lies in C, we see that PcA x (l)-Pc = ^ 2 Pc-> which implies 

a = A*(l)=p. (72) 

Stringing together the left inequality in Eq. ([7(]) and the right inequality in Eq. ( |7l| ) gives 

S(A P ) > S(p') - £ KS(pa) = S(p') - £ X a S(a a ) > S(p') - S(p) . (73) 

a a 

Condition 2 dictates equality all the way across here: equality on the left implies that the density 
operators p a are orthogonal, and equality on the right implies that the density operators a a are all 
equal to p. Thus we have that A p has the form fl55|), with orthogonal density operators p a . Since 
p has support on the entirety of C, we can put Ac in the reversible form (|53|), thus completing the 
proof. 



4.2 Bounds on entropy exchange and entanglement fidelity 

Thus far we have dealt with exact exact error correction. Both more difficult and more important is 
approximate error correction, for which we need a measure of the fidelity of a reversal. Schumacher 
[16] has introduced a suitable measure, 

F e {p,A)^^^\ P RQ '\^^) , (74) 

called the entanglement fidelity] it measures the fidelity with which A preserves the entanglement 



of the primary system with the reference system. Using the superoperator formalism of Sec. 3.2, 
we can write the entanglement fidelity as a superoperator matrix element: 

F e (p,A) = (^p\A p \y/p) = (p\A\p) . (75) 

Since A p is positive and has unit trace relative to the left-right action, F e (p, A) = 1 if and only if 
yfp is a (normalized) eigenoperator of A p with eigenvalue 1, i.e., A p = yfp®\fp~- In this case A p is 
a pure completely positive superoperator; both the von Neumann entropy (or entropy exchange), 
S(A P ) = 0, and the quadratic entropy, Tr(^) = 1, faithfully report this purity. To analyze 
approximate reversal, we need relations between these two measures of purity and the entanglement 
fidelity. Obtaining such relations is the task of this subsection. 

We look first at ways to use the purity measures to bound the entanglement fidelity away from 
1. Let the operators rj a , a = 1, . . . , D 2 , be an orthonormal operator basis chosen so that rj\ = ^/p, 
which means that F e = (r]i\A p \r]i). Then we have 

D 2 D 2 

S(A P ) < -^2(r] a \Ap\r] a )log(rj a \Ap\r] a ) = -F e logF e - ^ {iq a \Ap\r} a ) \og(r] a \A p \r} a ) . (76) 
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The inequality here is the standard result that the von Neumann entropy is never greater than the 
entropy calculated in any orthonormal basis [25, 26]. Renormalizing the distribution remaining in 
the sum, we can write 

S(A P ) < h(F e ) + (1 - F e ) {- f M^bM log (vMM ^j f (77) 

where h{x) = — xlogx — (1 — x) log(l — x) is the binary entropy. Using the fact that the entropy 
within the large parentheses is bounded above by log(D 2 — 1), we get Schumacher's [16] quantum 
Fano inequality, 

S(Ap) < h(F e ) + (1 - F e ) \og(D 2 - 1) , (78) 

which bounds the entanglement fidelity away from one. Schumacher obtained the quantum Fano 
inequality by applying the same reasoning to p R ® . 

We can do the same thing with the quadratic entropy. The analogue of Eqs. (^) and (|77|) is 



Tr(^) = £|MA>M| 2 > Efa»WK) a = F " + (1 " E 

aS a=l a=2 \ 



(79) 



Since the remaining sum is bounded below by (D 2 — l) -1 , we get a quadratic quantum Fano 
inequality, 

Tr (^) > F " + ' (80) 

which like the entropy version, bounds the entanglement fidelity away from one. 

We now turn to the opposite task: using the purity measures to place lower bounds on the 
entanglement fidelity. That this might work is suggested by the case S(A P ) = or, equivalently, 
Tr(v4p) = 1. Then A p = r (g> -r is pure, the normalized eigenoperator r having eigenvalue 1. If A is 
trace preserving, we have that p = Ap (1) = rV. Then, by the polar-decomposition theorem, there 
exists a unitary operator U such that t = U^fp. Defining a new operation by 

A = U ] ®U oA, (81) 

we have that A'„ = yfp®yfp, so that F e (p,A') = 1. What we have shown is that A is within a 
unitary of an operation that has unity entanglement fidelity for input density operator p. 

We now mimic this construction when A p is not pure. In doing so, we assume that A is 
trace preserving. Suppose that A p has the eigenvalue decomposition (|6"6|). Let Ai be the largest 
eigenvalue, and let 



n = U^t\t x = Uyfd x , (82) 
as in Eq. (p7|). Define the new operation (|8l|), with the result that 



A' p = Ai v /oT«>v / ar+ E X ^r a t\U . (83) 



Now we use the fact that A is trace preserving to write 



p = A* (1) = Aiai + Kr ] a T a , (84) 
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which implies that p > Aioi. Because the operator square-root function is an operator-monotone 
function [[^j, this operator inequality remains true upon taking the square root of both sides: 



Vp>V\iV^- (85) 
The operator-monotone property of the square root is is proved in the Appendix. By writing o~i in 



terms of its eigendecomposition, we see that the operator inequality (|85|) implies that 

tr(VPv^D > v^tr^i) = VAT • (86) 

Now we notice that 

F e (p,A') = (VpI^IVp) > \i{tr{^p^~i)f > A? • (87) 

This is our key result. It says that if the largest eigenvalue of A p is close to 1, then A can be 
corrected by a unitary so that the entanglement fidelity is close to 1. 

We can translate Eq. (^7j) into weaker bounds that involve the purity measures. For the entropy 
exchange we have that 

S(A P ) = - A « log A « > " E A ° lo S A i = " lo S A i • (88) 

a a 

Using Eq. (p7\), we obtain a quantum anti-Fano inequality 

F e (p,A')>e W (-2S(A p )) . (89) 
For the quadratic entropy we can write 

^p) = E A ^E A - A i = A i- ( 9 °) 



Again using Eq. (|87|), we obtain a quadratic quantum anti-Fano inequality 



F e (p,A')>(Tv(Al)) 2 . (91) 



Both of the anti-Fano inequalities place lower bounds on the entanglement fidelity as we had hoped. 
It should be remembered that they apply only to trace-preserving operations. 

Appendix 

Let A and B be positive operators such that A > B. Define the Hermitian operator 

F = VA-VB = J2Fk\fk)(fk\, (92) 

k 

where the sum is the eigenvalue decomposition of F. The eigenvalues of F satisfy 

F k = {f k \VA\f k ) - (f k \VB\f k ) > -(f k \VB\f k ) . (93) 
Now define the Hermitian operator 

G = A-B = VBF + F^B + F 2 (94) 

Since G > 0, we have 

< </*|G|/*> = 2F k (f k \VB \f k ) + Fl . (95) 
Equations (j93| ) and ( p5|) together imply that F k > 0, which means that VA > \f~B. 
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Table 1: How the 27 single-qubit error operators affect an arbitrary superposition of the logical 
zero \0)l = |TTT) and the logical one 1 1) x, = in Shor's nine-bit code. The first line shows the 
superposition state with no error. 
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